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Abstract. We present the most suitable data sets available in the International Ultraviolet Explorer (IU E) archive 
for the study of time-dependent stellar winds in early B supergiants. The UV line profile variability in 11 B0 to B3 
stars is analysed, compared and discussed, based on 16 separate data sets comprising over 600 homogeneously re- 
duced high-resolution spectrograms. The targets include ’normal 1 stars with moderate rotation rates and examples 
of rapid rotators. A gallery of grey-scale images (dynamic spectra) is presented, which demonstrates the richness 
and range of wind variability and highlights different structures in the winds of these stars. This work emphasises 
the suitability of B supergiants tor wind studies, under-pinned by the fact that they exhibit unsaturated wind 
lines for a wide range of ionization. 

The wind activity of B supergiants is substantial and has highly varied characteristics. The variability evident in 
individual stars is classified and described in terms of discrete absorption components, spontaneous absorption, 
bowed structures, recurrence, and ionization variability and stratification. Similar structures can occur in stars of 
different fundamental parameters, but also different structures may occur in the same star at a given epoch. We 
discuss the physical phenomena that may be associated with the spectral signatures, and highlight the challenges 
that these phenomena present to theoretical studies of time-dependent outflow's in massive stairs. 

In addition, SEI line-synthesis modelling of the UV wind lines is used to provide further information about the 
state of the winds in our program stars. Typically the range, implied by the line profile variability, in the product 
of mass-loss rate and ion fraction (M q t ) is a factor of ~ L5. w r hen integrated between 0.2 and 0.9 it can 
however be several times larger over localised velocity regions. At a given effective temperature the mean relative 
ion ratios can differ by a factor of 5. The general excess in predicted (forward-scattered) emission in the low 
velocity regime is discussed in terms of structured outflows. Mean ion fractions are estimated over the B0 to Bl 
spectral classes, and trends in the ionic ratios as a function of wind velocity are described. The low' values obtained 
for the ion fractions of UV resonance lines may reflect the role of clumping in the wind. 
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1. Introduction 

The radiation-pressure-driven winds of O and early-B type 
stars are variable on characteristic time-scales of hours and 
days, as a result of perturbations in the outflow due to 
small- and large-scale structure, respectively. Observable 
diagnostics of the former include extended black absorp- 
tion troughs in saturated UV P Cygni profiles (Lucy 
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1983). X-rav emission (Feldmeier et al. 1997, Kahn et al. 
2001, Cassinelli et al. 2001), and non-thermal radio emis- 
sion (Bieging et al. 1989, Scuderi et al. 1998). The small- 
scale wind structures responsible for these signatures are 
believed to result from a strong instability in the line- 
driven winds, which operates on linear scales smaller than 
the Sobolev length. This instability leads to multiply non- 
monotonic velocity fields in the wind and localised density 
enhancements which are bounded by pairs of forward and 
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reverse shocks (e.g. Ovvocki et aL 1988: Feldmeier 199b) . 
Large-scale structure in hot star winds is primarily stud- 
ied via systematic variability of UV resonance and opti- 
cal emission lines. These changes diagnose organised, spa- 
tially extended perturbations in the wind, which may be 
induced by magnetic or pulsational processes arising at the 
stellar surface. We document in this paper the nature of 
large-scale wind structure in 11 early-type B supergiant 
stars, based on high-resolution International Ultraviolet 
Explorer (IU E) time-series spectroscopy. 

One of the major discoveries made with II : E was the 
systematic nature of variability in the P Uvgiii-type line 
profiles of OB stars. The key developmental sequence was 
defined in a number of time-series UV investigations (e.g. 
Prinja et al. 1987, Prinja & Howarth 1988. Henrichs et al. 
1988, Kaper 1993, Kaper et. al. 1996). A relatively broad, 
shallow absorption feature appears at < 0.3 of the terminal 
velocity (Cog)i which then evolves blue- wards over several 
hours and perhaps even days, while narrowing in velocity 
width. The velocity at which these 'discrete absorption 
components 1 (DACs) are first detected is roughly constant 
for a given star, but may vary from star-to-star. 

Subsequent IU E studies established a connection be- 
tween stellar projected rotation velocity ( v e sin(z)) and 
the rate of recurrence and acceleration of DACs (Prinja 
1988; Kaper et al. 1996). More recent IU E monitoring 
campaigns extending over a few weeks have shown that 
the stellar wind lines vary continuously, implying that OB 
star winds are structured over a range of spatial and tem- 
poral scales. They have also shown that variations in the 
profiles modulate on rather long time-scales (/.<?. several 
days). Furthermore, the variability evident in £ Per (07.5 
111(f); Henrichs et al. 1998) and £ Pup (04 If; Howarth et 
al. 1995) appears to modulate on the stellar rotation rate. 
Rotation is plausibly also the key factor in the repetitive 
wind activity seen in HD 64760 (BO. 5 lb; Fullerton et al. 
1997), HD 93843 (05 III; Prinja et al. 1998) and HD 91969 
(BO la; Massa et al. 1998). Rotationally modulated wind 
fluctuations are also apparent in late B and early A super- 
giants (e.g. Kaufer et al. 1996), which have much larger 
radii and longer rotation periods. 

The fundamental origin of wind variability on large 
spatial scales in hot stars remains uncertain. The causal 
link between inhomogeneities at the stellar surface and 
large-scale wind structure has not been directly estab- 
lished observationally or theoretically. A currently pop- 
ular idea is that the large-scale variations discussed above 
(and some of those discussed in this paper) are spawned 
by photospheric irregularities which cause the wind from 
different longitudinal sectors on the surface to emerge with 
different densities and/or velocities. The surface structure 
may, for example, be due to the action of non-radial pulsa- 
tions or the effects of ordered magnetic fields. Once differ- 
ent adjacent wind streams are initiated, they then meet to 
form co-rotation interaction regions (CIRs). The rotating 
CIRs inevitably evolve into spiral-shaped shock structures 
(Cranmer & Owocki 1996). In this model, variations in 


the wind absorption result from material flowing through 
a semi-permanent, co-rot ating wind structure. 

1.1. Winds of early B-type supergiants 

The stellar winds of early B supergiants (Bis) are par- 
ticularly suitable for a study of wind dynamics and the 
effects of large-scale wind structure since, (i) their UV 
wind lines are often well-developed but unsaturated, (ii) 
these diagnostic lines cover a w'ide range of ionization, 
close to the dominant stage of the wind (see Table 2), 
(iii) it is generally easier to follow' the progression of 
large-scale structure in the UV lines of Bis compared 
to O-type stars, (iv) the winds have lower terminal ve- 
locities, so individual structures evolve more slowly and 
the important Siiv AA1400 resonance line doublets are 
mostly decoupled, (v) their strong UV photospheric lines 
are excellent effective temperature-luminosity diagnostics, 
and (vi) the strongest UV photospheric lines sample the 
wind/photosphere interface 1 , often displaying wind be- 
haviour. 

Motivated by these attractions, we describe in this 
paper the time-dependent wind characteristics of eleven 
BO to B3 supergiants, based on 16 separate IUE high- 
resolution time-series data sets, involving the inspection of 
638 spectrograms in total. This work extends our knowl- 
edge of the characteristics of variable UV resonance lines 
in line-driven winds to cooler and less luminous stars. The 
data are used to demonstrate the presence of a range of 
empirical wind properties and implied varied perturba- 
tions in the outflows. We report here on the incidence 
and repetition of DACs, evidence for rotational modula- 
tion, ionization-state changes, shock formation, and spon- 
taneous wind enhancements. The diverse wind variability 
signatures in early Bis are compared and contrasted to 
those previously noted in luminous O-type stars. The later 
sections of this paper present SEI line-synthesis modelling, 
which is used to explore further the ionization mixture and 
structured nature of the winds. 

2. The IUE time-series data sets 

Our objective was to carry out a systematic survey of the 
temporal nature of wind activity in early B supergiants. 
We therefore restricted our selection to time-series of high- 
resolution (A/AA - 10 4 ) IUE SWP (AA1150 to 1900A) 
data sets. A list of the program stars is given in Table 1. 
They include luminosity classes la and lb for a spectral 
range BO to B3, and have a wide range of projected ro- 
tation velocities. The corresponding 16 data sets analysed 
are also outlined in Table 1, and include 4 cases where 
short-term wind variability may be compared over two or 
more epochs separated by a few years. The majority of 
stars in our sample are 'normal 1 , slow to moderate rota- 
tors, that are representative of their spectral sub-class. 
Two main exceptions are HD 64760 and HD 157246 (7 
Ara) which were selected to examine the effects of rapid 
rotation on wind structure. 
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Table 1. Program stars 


HD 

Sp. Type Log( L. i L 

,) Tnff [kK\ r r . sin(/) 

1 km s' 1 ) 

Protfinax) 

(days) 

Date of IUE 
t ime series 

Duration 

(days) 

Sampling 

(hours) 

37128 

BO la N“ 

5.78 

28.6 

91 

17.7 

1980 

0.7 

1 

91969 

BO la 

5.78 

28.6 

83 

14.4 

1996 

30 

9 

164402 

B0 lb N" 

5.38 

28.6 

77 

13.3 

1986/93/95 

6/16/16 

3-24/24/21 

204172 

B0 lb N" 

5.38 

28.0 

87 

11.6 

1990/93 

2/16 

3/24 

64760 

BO. 5 lb 

5.10 

23.1 

216 

5.2 

1993/95 

6/15 

3/3 

167756 

BO. 5 lb 

5.50 

23. 1 

79 

14.1 

1986 

5.5 

12 

150168 

Bl lab-lb 

5.38 

20.3 

121 

13.3 

1993 

6 

3 

47240 

Bl lb 

4.90 

20.3 

103 

1 1.3 

1993 

16 

24 

157246 

Bl lb 

4.90 

20.3 

245 

4.8 

1993/95 

6/16 

3/24 

96248 

BCl la 

5.12 

19.4 

83 

24.3 

1996 

30 

9 

53138 

B3 la 

5.24 

16.3 

58 

45.8 

1987 

13 

12-36 

and Log(L,/L;,) from Hum 

phrevs 

A Me: Elroy 

f 1984); tv sin( t) 

from Howarth et al. (1997) 

Spectral types from VValbom 

1972, 1976), 

Garrison et al. ( 

1977), 

Hiltner et al 

( 1969), and Morgan et al. 

(1955). N refers to stars 

whose spectra are 


nitrogen deficient. 


The data were either extracted from the archives main- 
tained at the Goddard Space Flight Center (I SA) and 
the Rutherford Appleton Laboratory (UK), or they were 
originally part of our observing campaigns. All the spec- 
tra were extracted in a uniform manner from the two- 
dimensional, photometrically linearised images (PI) using 
the UK Stariink package IUEDR (Giddings et al. 1995). 
One of the advantages of the IV EDR extractions is the 
spectra are highly reliable at the shortest wavelengths, 
thus providing high data quality for the important lines at 
Si ill A1206 and N v AA1240. The wavelength scale calibra- 
tion was improved by centering the spectra on selected in- 
terstellar lines (in an echelle order-dependent correction). 
The individual spectra were binned onto a wavelength grid 
with a regular sampling of 0.1 A. Typically the signal-to- 
noise in the continuum is ~ 20. Individual wind-formed 
line profiles were normalised to a local continuum by fit- 
ting a straight line through several adjacent line-free* win- 
dows. 

The mean IV E SWP spectrum for each program star 
is shown in Fig. 1. The principal diagnostic wind and 
photospheric lines available in these spectra are listed in 
Table 2 (Massa 1989). The data sets for the following stars 
have been previously reported in the literature: HD 64760 
(Massa et al. 1995b; Prinja et al. 1995: Fullerton et al. 
1997;Howarth et al. 1998), HD 157249 (Prinja et al. 1997), 
and HD 91969 (Massa et al. 1998). These data are pre- 
sented in this paper to place them into wider comparative 
context. 

3. Time-dependent stellar wind phenomena 

The characteristics of variable UV lines in our pro- 
gram stars are presented in this section. Fluctuations are 
present in the stellar wind lines of every B supergiant in 
our sample, and in every time-series data set examined 
here. The variations are primarily described by a gallery 
of grey-scale image representations (dynamic spectra) of 
the time-series data (Figs 2 to 17). This image format 
is an effective method for identifying systematic patterns 


in variable line profiles, and for comparing the evolution 
of these features between different spectral lines. In cases 
where the Doppler shifted absorption profiles are unsatu- 
rated (i.e. the majority in our sample), quotient spectra 
are displayed after being normalised by the mean pro- 
file for the time-series. For saturated or ‘black' absorption 
troughs in the Civ lines of HDs 37128, 91969, 157246 and 
96248, we display residuals, subtracted from the mean pro- 
file. In all cases darker shades in the images indicate re- 
gions of enhanced absorption (reduced flux) with respect 
to the mean profile. 

3.1. List of observed stellar wind phenomena 

Several empirical phenomena can be identified from the 
images in Figs 2 to 17. They are defined below, prior to 
discussion of their occurrence in individual stars: 

Discrete absorption components (DACs) - these are the 
features most commonly associated with variable UV res- 
onance line profiles in luminous OB stars. They usually 
appear as localised (in velocity) absorption enhancements 
that migrate from (at least) intermediate velocities to the 
blue-ward profile edge, with a corresponding decrease in 
velocity width. Typically the accelerations of the DACs 
are < 50% of values for the underlying wind, as deter- 
mined from canonical ‘/? ~ V velocity laws. 

Bowed structures - these are optical depth changes that 
appear suddenly over a wide velocity range, and which 
develop simultaneously towards lower and higher blue- 
shifted velocities (see e.g. Fullerton et al. 1997). 

Other features - we use this term to identify either less 
organised ‘bursts' of absorption that usually last a short 
time only (hours), or more substantial absorptions that 
appear suddenly over a large velocity range. These fea- 
tures do not show any clear evolution in their velocity or 
width as a function of time. 

Near-zero velocity changes - cases where line profile vari- 
ability can be traced to less than 0.05 of the terminal 
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Fig. 1. Mean normalised spectra for the IUE time-series of the 11 program stars. The important spectral lines are identified in 
Table 2. (An intensity shift of 1.5 has been used in this display between successive spectra.) 

velocity. These fluctuations are due to deep-seated distur- repetitive structures. The repetition is normally due to 
bances close to the stellar surface. the occurrence of DACs or bowed structures. 


Recurrence - This is simply a descriptive term relating Ionization variation - implies that the ionization state of 
to a specific (extended) data set, to draw attention to a specific wind feature changes with time. It may for ex- 
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Table 2. Key IE E line diagnostics foi early-B supergiants winds 


Ion 

Ao (A) 

I.P. (ev) 

E.P. (c-V) 

region 

(’ll 

1414.53. 1335.71 

2 1.4 

0.00. 0.01 

wind 

C iv 

15 10.20. 1550.77 

64.5 

0.0, 0.0 

wind 

Nv 

123X.X2. 1242.X0 

97.9 

0.0. 0.0 

wind 

Al ill 

1*54.72. 1X62.79 

28.5 

0.0, 0.0 

wind 

Si ill 

1206.50 

33.5 

0.0 

wind 

Si in 

1296.73 

33.5 

6.55 

wind / photosphere 

Si in 

1312.24 

33.5 

10.28 

photosphere 

Si in 

1417.24 

33.5 

10.23 

photosphere 

Si iv 

1393.76. 1402.77 

45.1 

0.0. 0.0 

wind 


ample become more enhanced in absorption in some lines, 
while exhibiting an increased flux in other spectral lines. 
Ionization stratification - cases where one of the phenom- 
ena identified above are seen at different velocities in dif- 
ferent ions. This is usually an optical dept h effect between 
low and high ionization lines. 

Results on the occurrence of these variable wind phe- 
nomena is our program stars are summarised in Fable 3. 
The winds of early Bis exhibit diverse patterns of stel- 
lar wind variability. Similar phenomena can occur in stars 
that vary in fundamental parameters, and different wind 
structures can co-exist in the wind of a given star at a par- 
ticular epoch. Some notes on the temporal wind behaviour 
in the individual stars are presented in the remainder of 
this section. 

3.2 . Notes on individual stars 

HD 37128 (e Ori, BO la N~; Fig. 2). A well-studied, mod- 
erately nitrogen weak star (Walbo n 1976), part, of the 
Ori OBI association, c Ori has been cited as a single-lined 
spectroscopic binary (e.g. Frost 1909; Morrell k Levato 
1991), however its binarity remains uncertain and, for ex- 
ample, Jarad et al. (1989) found it to be a radial velocity 
constant. Variability in the hydrogen lines has been re- 
ported by Gry et al. (1984) and Ebbets (1982). Previous 
suggestions that it is a variable thermal radio source are 
discounted by Blomme et al. (2002) 

UV variability in the short intensive IV E time-series is 
shown in Fig. 2. VVe note the progression of a DAC feature 

between 900 to —1200 km s' 1 ewer — 14 hours. High 

velocity narrow absorption components are also present 
towards the short-ward edge. Strong photospheric lines 
(e.g. Cm 1247 and Aim 1855) reveal radial velocity shifts 
of — 20 km s' 1 . 

HD 91969 (B0 la; Fig. 3). A member of the (Jar OBI 
association (e.g. Humphreys 1978). It is one of the bright- 
est members of the NGC 3293 open cluster in the Carina 
region (e.g. Feinstein k Marracu 1980). The IU E time- 
series data studied here were secured during our "MEGA 
IF campaign in 1996, and some of our findings have been 
summarised by Massa et al. (1998). This star is an excel- 
lent example of the variety of diagnostic UV lines available 
in early Bis. Massa et al. (1998) reported a 7.8-davs period 
for the repetition of DAC features propagating blue-ward 


toward the profile edge. These features are seen at very 
low velocities in Aim, and their progression ‘connects’ 
with corresponding events seen at higher velocities in Si ill. 
Si IV and Nv. There is however evidence for ion stratifi- 
cation between activity in Aim and Si IV. The temporal 
coherency is strongest in the low ion species which form 
in the deeper regions of the wind. The periodic line profile 
changes are seen down to very low velocities (less than 
0.05 I'coh together with variations in photospheric lines 
due to Sim, Fem, and Feiv. The latter may be tracers 
of localised surface inhomogeneities responsible for wind 
activity. 

HD 164402 (B0 lb N“; Figs. 4 to 6). A member of the Sgr 
OBI association (Humphreys 1978), this star has one of 
the lowest projected rotation velocities in our sample. It 
is moderately nitrogen deficient, with no indication of bi- 
nary companionship (e.g. Levato et al. 1988). HD 164402 
provides an opportunity to study systematic wind activ- 
ity over three separate epochs. The 1986 data set (Fig. 
4) has a poorer temporal sampling, but reveals a distinc- 
tive DAC progressing from — -1200 km s' 1 to ^ -1600 
km s _1 over 6 days. This feature is accompanied by more 
sporadic events, which are well spread in velocity, but un- 
resolved temporally (e.g. at T ~ 3.5, 4,2 and 5 days, where 
T is time relative to the first observation in the series). The 
1993 time-series (Fig. 5) spans almost 16 days and reveals 
a remarkably slowly migrating DAC, with a mean acceler- 
ation of ~ 3 x 10~ 4 km s’ 2 . Other features with different 
characteristics and more rapid evolution are also noted at 
T ~ 6 and 12 days. Drastic variations are present in the 
1995 data set (Fig. 6) around T ^ 14 days. The high- 
est velocity absorption (at — — 1750 km s" 1 ) suddenly 
decreases and is accompanied by lower velocity DAC fea- 
tures at ^ —1300 km s’ 1 and —1000 km s' 1 . The sudden 
reversal in optical depth across the lines may represent 
the signature of a large scale shock. 

HD 204172 (69 Cyg; B0 lb N“: Figs. 7,8). A moder- 
ately nitrogen weak star (Walborn 1976), belonging to 
the Cyg OB2 association. Lennon et al. (1992) highlight 
several discrepancies between the B0 lb spectral type as- 
signment and the appearance of the optical spectral line 
profiles, including H lines being too narrow and Si IV too 
strong. They suggest a revision to B0. 2 la. The UV reso- 
nance lines of 69 Cyg are generally stronger than those of 
HD 164402, thus favouring a slightly more luminous spec- 
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Table 3. Observed stellar wind phenomena 


HI) 

DAG 

Bowed 

structure 

Other 

feature 

near- zero 
var. 

ion 
st rat. 

ion 

var. 

recurrence 

37128 

V 


- 

- 

- 

- 

- 

91969 

V 


y 

y/ 

y/ 

y 

y 

164 '10 2 

V 


V 

O 

- 

- 

- 

204172 

■ > 


V 7 

2 

... 

- 

- 

64760 



v/ 

sj 

yj 

y 

v/ 

167756 

V 

- 

V 

yj 

yj 


- 

150168 

V 

v/ 

y/ 

y/ 

yj 

y 

v' 

47240 

v/ 

- 

y/ 

si 

yj 

- 

- 

157246 

y/ 

- 

yj 

- 

yj 

y 

- 

96248 

v/ 

- 

v/ 

- 

yj 

- 

- 

53138 

sj 


y/ 

- 

- 

- 

- 


tral type. The densely sampled 1990 time series (Fig. 7) 
reveals a gradual increase in profile absorption strength 
over — 2.2. days. The fluctuations extend over a broad ve- 
locity range that cover almost the entire Siiv absorption 
trough. Corresponding changes are seen at much lower ve- 
locity (extending close to zero velocity) in Cn A1335 and 
Aim A 1855. The velocity differences between low and high 
ions are principally an optical depth effect. There are also 
traces of DAC features at ^ —1700 km s -1 (e.g. T 0 to 
0.5 days) in Si IV, C iv and N v. Some changes are appar- 
ent in the strong Cm A1247 "photospheric' line. The 1993 
data set (Fig. 8 ) is more extended and covers almost 16 
days. Two enhanced absorptions are present over a sub- 
stantial velocity range. These may be similar events to the 
phenomena seen with greater time resolution in 1990 (Fig. 
7). There is a tentative indication that the feature recurs 
over ~ 12 days, which is close to the estimated maximum 
rotation period of the star (Table 1). Some variations are 
also noted in the deep-seated C III A 1247 and Si III A 1300 
triplets. 

HD 64760 (B0. 5 lb; Figs. 9, 10). This is one of two rapid 
rotators in our sample. The data sets for the two epochs 
shown here were discussed in detail by Massa et al. (1995b) 
and Fullerton et al. (1997). HD 64760 exhibits a variety of 
spectral features that indicate co-existing wind structure. 
In fact, every category of stellar wind phenomena listed in 
Table 3 is present in this star! The "MEGA F campaign 
during 1995 (Fig. 10) in particular reveals a striking 1.2- 
day modulation, where individual features cover a large 
velocity range, and have a bow-shaped morphology such 
that they extend blue-wards and red-wards at the same 
time. Ionization variability is indicated by the fact that the 
N v data appear to have nearly twice the modulation fre- 
quency of Si IV (Fullerton et al., 1997). Slower migrating, 
high velocity DACs are also present, which persist over 
several days. The recurrence time-scale for the DACs is 
not constrained. An exceptionally strong absorption event 
is present in high-ion species in the (shorter) time series 
data for 1993 (Fig. 9, around T ~ 2 days). It is accom- 
panied by variations essentially down to zero velocity in 
e.g. C II AA1335 and Aim A1855. plus disturbances in the 
‘photospheric 1 Si m A1300 triplets. The onset of these sub- 


stantial changes is also marked by changes in the relative 
ionization mixture, such that the overall ionization of the 
wind increases during this event (Massa et al. 1995b). It 
is possible that this strong event, and the overall enhance- 
ment seen at T ~ 4 days (Fig. 9), are similar to the bowed 
structures apparent in the more extensive MEGA I data 
set (Fig. 10). 

HD 167756 (B0. 5 lb; Fig. 11). A distant halo star, located 
about 4 Kpc away at a Galactic altitude of —0.85 Kpc 
(e.g. Savage & Massa 1985). Ionization shifts are appar- 
ent in the wind-formed lines, with the high velocity 

— 1500 km s" 1 ) DAC seen in Siiv (T ~ 0 to 3 days) only 
being apparent at less than —500 km s -1 in Aim. The ab- 
sorption enhancement seen in Siiv, N v, and Civ towards 
the end of the observing run (T > 4 days) may also be 
accompanied by variations close to zero velocity in Sim 
A1207 and Aim A1855. 

HD 150168 (B1 Iab-Ib; Fig. 12). A member of the Ara 
OBI association (Humphreys 1978). The UV time-series 
reveals a radial velocity motion in strong UV photospheric 
lines (e.g. Cm A1247) due to binary motion. Several in- 
teresting wind effects are apparent in this star. Near the 
start of the time-series there is a maximum in absorption 
in N v which corresponds to a minimum in Siiv (i.e. ion- 
ization variability). Close to T ~ 4 days a bowed structure 
is present, with the absorption minimum in N v lagging 
the Siiv minimum (i.e. ion stratification). A much weaker 
bowed feature may be present at T — 2 days. Multiple 
high velocity DACs (> —1500 km s -1 ) are seen in Siiv, 
Civ and N v. 

HD 47240 (B1 lb; Fig. 13). A member of the Mon OB2 
association. Binary radial velocity motion is seen in strong 
UV photospheric lines (e.g. Cm A1247).The double DAC 
feature seen at the start of the time-series at ~ —700 km 
s _1 in Siiv is present at a lower velocity (< —500 km s~ l ) 
in C II and Al III. This enhancement is seen at almost zero 
velocity in Civ, however there is very little line flux in 
Si III A1207. The variations in Siiv between — —500 and 

— 1000 km s -1 may be repetitive over a 3 to 5 day time 
scale. 

HD 157246 (7 Ara, B1 lb; Figs. 14,15). This star has 
the highest projected rotation velocity in our sample of 
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stars ( Fable 1 ). The 1993 UV tirne-sei ies was discussed by 
Prinja et al. (1997: see Fig. 14). T he data reveal a com- 
plex pattern of co-existing discrete features at. low veloci- 
ties (red- ward of —750 km s -1 ) and high velocities (blue- 
ward of —1500 km s _1 ). These region s are ‘divided’ by 
the appearance of a very sharp (in velocity width) DAC- 
like feature. The location of this narrow feature is shifted 
by — 400 km s " 1 between Si IV and Aim (or Cm). The 
two-component wind properties, and overall profile mor- 
phologies, suggest, evidence for an equ itorially compressed 
wind, with structure in the polar components (Prinja et al. 
1997). The narrow DAC-like feature is not seen in the 1995 
IU E time-series (Fig. 15), and the high speed wind be- 
tween — 1000 and —1500 km s“ l is also no longer present. 
The line profiles in the 1995 data are more representa- 
tive of the normal state of 7 Ara (e.g. when compared to 
isolated archival spectra taken at earlier epochs). Some 
ionization state variability is observed between Si IV and 
N v (see e.g. the fluctuations between T — 2 to 8 days in 
Fig. 15). In both data sets variations are present down to 
near-zero velocities in the low ions. 

HD 962^8 (BC1 la; Fig. 16). Carbon enhancement has 
been reported in this star by Dufton (1972) and Walborn 
(1976). HD 96248 is a radial velocity (onstant (e.g. Levato 
et al. 1988). The ~ 30 day time-series shown in Fig. 16 is 
from the 'MEGA IT Campaign of 1996. The Si IV and C IV 
profiles are nearly saturated, and the line profile variability 
is confined to the blue edges (i.e. > —750 to —1000 km 
s" 1 ). A gradual increase in absorption in the blue edge 
occurs from the start of the time-series to almost the end 
of the run, only weakening slightly after F ~~ 24 days. 
These variations are mimicked in the low ions (Aim and 
Cll), but at slightly lower velocities. Additional variations 
due to DACs are also seen between ^ —200 and —600 km 
s -1 in Aim, but they do not evolve into the structure seen 
at high velocities. The blue edge in Cll A1335 is not as 
shallow or extended as in the other lines. 

HD 53138 (o 2 CMa, B3 la; Fig. 17). A member of the Coll 
121 association. HD 53138 exhibited very complex and 
highly variable characteristics in the Ha profile survey 
of Ebbets (1982), with variable absorption and emission 
components. Low velocity (red- ward of —500 km s” 1 ) UV 
variability in present in the ~ 13 day time series shown 
in Fig. 17. A DAC feature is preseni between T ~ 0 to 
4 days, which is essentially confined to high velocities. 
Various sporadic changes are apparent, later in the run. 
Some degree of ion stratification is present, such that Cll 
variations during the first 4 days appear at lower velocities 
than in Al III. 

4. Comparisons 

It is clear from Figs. 2 to 17 (and Table 2) that the stellar 
wind activity in early B supergiants as revealed by the 
UV resonance lines, is substantial ai d has highly varied 
characteristics. 

It is interesting to compare the Bl wind activity pre- 
sented here with that of 10 O-lype stars described by 


Kaper et al. (1996), also based on high-resolution IV E 
spectroscopy. The O stars generally do not exhibit such 
a variety of temporal phenomena, and in particular show 
limited evidence for ionization related effects. (The com- 
parisons are restricted, however, since Si IV is often the 
only reliable indicator of wind activity in O-type stars, 
as N v and Civ are saturated and the sub-ordinate N IV 
1718 is only weakly present.) The fundamental character- 
istic patterns of variability in the O star data are in the 
form of blue-ward migrating DACs, though their proper- 
ties may differ a lot between e.g. acceleration rates, mul- 
tiplicity, and line-of-sight velocity dispersions. The coher- 
ence of systematic structure is generally greater in the 
UV lines of Bis compared to those of O stars. In terms of 
clearly defined - obviously modulated - wand structure, 
the B supergiant, examples of HD 64760 and HD 91969 
described above, are best matched by the 07.5 III stars £ 
Per (HD 24912) and 68 Cygni (HD 203064). 

In order to quantify the extent of variability as a func- 
tion of velocity, we computed the temporal variance spec- 
trum (TVS, see Fullerton et al. 1996) for the Si IV line 
profile for all our program stars: 


(TVS),- - 


o 




jV 

E 


( Sij - Si \ 2 

\ a jc\/Sij J 


( 1 ) 


where Sij is the normalized intensity of the ith pixel 
in the j th spectrum, Si is the weighted mean of the 
normalized intensity, <jj c is the inverse of S/N of spec- 
trum j measured in an adjacent continuum band, and 
^ r | _ ^ ^ ^ 1 1 ^ 

a o — kv 2 ^j = i I • The results are shown as root 

mean square percentages (= TVS^ x 100) in Fig. 18. In 
all the stars with extended time-series data (i.e. spanning 
more that a few days), the absorption variability in Siiv is 
observed to extend down to less than 0.1 and to essen- 
tially rest velocity in a few cases (e.g. HDs 64760, 150168, 
47240 and 157246). This is a lower velocity than is gener- 
ally the case for O-type stars, and the result suggests that 
the winds of B stars are structured differently. Greater 
variability in the lower velocity portion of the P Cygni 
profiles in Bis indicates a lack of forward-scattered emis- 
sion compared to 0 stars. This could be due to a mul- 
tiply non-monotonic velocity law (which implies a higher 
degree of structure) or large-scale deviations from spheri- 
cal symmetry (which decrease the net amount of forward 
scattering) or clumping (less scattering in all directions). 

The maximum amplitude of the variations in Fig. 18 
mostly occurs between ~ 0.7 to 0.9 v <*,. HD 91969 and 
HD96248 exhibit the largest amplitude changes in Siiv. 
Though the variance diagnostics in Fig. 18 reveal no case 
of substantial variation in the emission components of the 
P Cygni profiles, weak systematic changes have been noted 
in the UV emission lobes of HD 64760 (e.g. Fullerton et 
al. 1997). 

Measurements were also made of the ratios of mini- 
mum flux to maximum flux as a function of velocity for 
selected lines in the time-series data of each star. The over- 
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Fig. 2. Grey-scale representation of variability in HD 37128. 
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Fig. 3. Grev-scale representation of variability in HD 91969. 
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y-scale representation of variability in HD 164402 (1986). 
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Fig. 5. Grey-scale representation of variability in HD 164402 (1993). 
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y-scale representation of variability in HD 164402 (1995). 
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Fig. 7. Grey-scale representation of variability in HD 204172 (1990). 
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-scale representation of variability in HD 204172 (1993). 
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Fig. 9. Grey-scale representation of variability in HD 64760 (1993). 
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Fig. 12. Grey-scale representation of variability in HD 150168. 
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Fig. 13. Grey-scale representation of variability in HD 47240. 
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Fig. 14. Grev-scale representation of variability in HD 157246 (1993). 
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Fig. 15. Grey-scale representation of variability in HD 157246 (1995). 
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Fig. 17. Grev-scale representation of variability in HD 53138. 
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Fig. 18. Temporal variance (RMS) for the Si IV spectral lines in our program stars. The velocity scale is with respect to the red 
component of the doublet, and the rest wavelength of the blue component is marked by an arrow. The horizontal dashed lines 
mark the level of significance for a probability of 95%. 

all trend in N V, Si IV and Civ is of a gradual decrease in to 0.68 (at 0.7 Cco), F m i n (v)/F max (u)[Si iv] ~ 0.80 (0.2 

Fmm ( c)/F max ( t' ) (i.e. greater relative absorption or opti- v^) to 0.33 (0.7 t^), and F min (t’)/F max (t;)[C iv] — 0.53 

cal depth) with increasing velocity. The mean ratios (over (0.2 Voo) to 0.05 (0.7 ty^). The most extreme variations in 

all stars) vary as: F mm (i;)/F max (tO[N ! v] - 0.81 (at 0.2 r N ) lines other than Si iv (Fig. 18) occur in HD 150168 and 
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HD 47240 in N v and C IV. The weak ?r A! Ill flux ratios do 
not exhibit a clear trend in velocity and mean values vary 
between ^ 0.75 to 0.9 - stronger absorption variations are 
seen in this line in HDs 47240, 157246 and 96248. 

4.1. Rotational modulation 

Aside from the data sets of HD 91969 (Fig. 3) and 
HD 64760 (1995; Fig. 10), the remaining WE time-series 
data of the program stars do not reveal unambiguous evi- 
dence for repeatable or cyclic modulations. In several cases 
this is likely because the WE observations are not ex- 
tended enough to cover stellar rotational time-scales. We 
nevertheless performed a Fourier period-finding analysis, 
using the CLEAN algorithm (Roberts et al. 1997) to itera- 
tively remove the window function from the Fourier spec- 
trum. We examined all the N V, Si IV and Civ time-series 
data sets for the program stars (where the line profiles had 
been previously normalised to the local continuum). The 
results indicate that (excluding HD 91969 and HD 64760) 
the only other cases of marginal significance where power 
is apparent across a reasonable fraction of the absorption 
trough (i.e. over few hundred km s“ 1 ) are (i) a modula- 
tion of ~ 4 days in the 1993 data ser, of HD 204172 (Fig. 
8), and (ii) fluctuations on ~ 1 day in the N v data set of 
HD 150168 (Fig. 12). 

Variability in the Siiv lines is compared more directly 
in the images shown in Fig. 19. where the stars are or- 
dered in terms of estimated maximum rotation period. 
There is perhaps some overall similarity between the mod- 
erate rotators HDs 47240, 204172, 164402 and 91969. Only 
HD 150168 shows bowed features that resemble the strik- 
ing patterns seen in HD 64760. There is also a resem- 
blance between the very slowly propagating, high veloc- 
ity DAC seen in HD 64760 (blue-wards of — —1200 km 
s” 1 ) and DACs seen towards the short-ward profile edge 
in HD 204172 and HD 164402. 

In O-type stars the UV resonance line profile variabil- 
ity is usually dominated by the blue-ward migrating, re- 
current, DACs. The apparent connection between modu- 
lated wind activity and stellar rotation in these stars is 
then based on measured trends in the recurrence and evo- 
lution time-scales of the DACs. We have shown in Sections 
3 and 4 that the line profile behaviour of Bis is more com- 
plex, such that the variance (e.g. Fig. 18) is generally not 
due to a single phenomenon such as DACs, but instead 
reflects the contributions of different variable structures. 
Furthermore, the empirical properties of these co-existing 
features do not all follow the same evolution time-scales. 
Examples of fairly localised (i.e. velocity widths — 0.1 to 
0.3 Vro ) blue-ward migrating features that resemble DACs 
in O-type stars are seen in our sample in HDs 91969, 
164402, 64760 (1995), 150168, 157246 (1993) and 96248. 
The measured mean accelerations for the DACs in the 
highest v e sin(i) stars, HDs 64760 and 157246, are — 90 
km s -1 d~ l and 140 km s” 1 d” 1 . respectively. The corre- 
sponding values for the slower rotator?, HDs 96248, 150168 


and 164402 are ~ 7 kin s -1 d” 1 . 35 km s -1 d _1 and 20 
km s _l cl” 1 . Therefore, aside from the other variability 
phenomena identified in Section 3.1, these results are con- 
sistent with the view that the time-scale of stellar wind 
activity in radiation-pressure-driven, luminous early-type 
stars is related to the rotation period of the underlying 
star. However, our study also highlights other broad ab- 
sorption features which can arise and evolve substantially 
over comparatively short time-scales (— 1 to 2 days) in 
seemingly moderate or slow rotators, such as HD 204172 
and HD 150168 (see e.g. Fig. 19). 

5. Profile fitting with the SEI code 

In this section we present the results from wind line pro- 
file modelling, to derive further details of the UV vari- 
ability characteristics and ionization states of early B su- 
pergiants. Specifically, our objective was to derive values 
for the product of mass-loss rate and ion fraction (M qi ), 
for the mean line profiles of the time-series and for cases 
of overall maximum and minimum observed absorption. 
We also examined trends in the relative ionization ratios. 
Our approach in this analysis is limited to the simplified 
method of employing steady-state line synthesis models 
to match observed profiles that are obviously affected by 
temporal wind structure and inhomogeneities in the out- 
flow. Nevertheless, some useful information can still be de- 
rived from the success and failures in matching line profiles 
due to different ions. 

We modelled the N V, C iv, Si IV and Si III (A1207) wind 
lines (and in a few cases also Aim and Cll) using the 
‘Sobolev with exact integration' (SEI) method described 
by Lamers et al. (1987). In addition, we adopted the modi- 
fied treatment of the wind optical depth law introduced by 
Massa et al. (1995b), such that the radial optical depth 
T 'rad(^) is treated in 10 independent velocity bins (each 
0.1 i>co wide). This approach provides greater flexibility in 
matching absorption optical depths in a wind that is af- 
fected by time-variable structure. We assume a standard 
parameterised law for the wind of the form 

w — Wq 4- (1 - w 0 )( 1 - l/x)* 3 , (2) 

where w = v/v ^ and x = r/R, where R is the stellar ra- 
dius. We select 0 = 1.0 and u?o = 0.01, and generally find 
that the model profiles do not differ grossly for values of 
d between 0.5 to 2.0. A range of 0.05 to 0.1 was adopted 
for Uturbi which is the small-scale velocity dispersion pa- 
rameter (i.e. the localised velocity dispersion is typically 
-lOOkms” 1 ). 

Since the UV photospheric spectrum of early Bis can 
be fairly strong, it is necessary to allow for their effects 
in the SEI profile fitting. The high-resolution WE pho- 
tospheric spectra of the following stars were adopted (by 
spectral type) as ‘standards’ and explicitly included as 
a lower boundary for the calculations: HD 36512 (B0 
V), HD 36960 (B0. 5 V), HD 31726 (B1 V), HD 36959 
(B1.5 V), HD 35468 (B2 III), HD 207330 (B2.5 III) and 
HD 209008 (B3 III). These 'standard* stars were selected 
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Fig. 19. Comparison of the Si IV variability characteristics of selected program stars; the estimated maximum rotation period 
is increasing from left to right. 


as low v e sin(i) targets, covering the range of spectral 
types of our program stars, and exhibiting very weak or no 
evidence for winds in their resonance lines. The input pho- 
tospheric spectra were convolved with a rotational profile 
to match the rotational broadening of the program stars, 
and generally the agreement between the corresponding 
photospheric absorption lines was very good. 

With the parameters in the velocity law and photo- 
spheric input fixed, the r rad (m) bins were the only free pa- 
rameters in the line fitting scheme. The procedure adopted 
was to fit the individual profile, bin by bin, progressing 
from the blue-ward profile edge to less negative veloci- 
ties (see Massa et al. 1995b). We modelled in this manner 
the mean Nv, Civ, SilV, and Si III profiles for the time- 
series data of each program star. In addition, to gauge 
the range in parameters due to profile variability, we also 
modelled for N v, Si IV, and Civ, the representative (over- 
all) maximum and minimum absorption spectra (the IUE 
SWP images used in these cases are listed in Table 4). 
The Sim A1207 lines are mostly too weak and their pho- 
tospheric absorption lines are less well matched by the 
standard stars, to determine variability characteristics by 
profile fitting. The mass-loss rate results are summarised 
in Table 4, where the M qi values are integrated over 0.2 


to 0.9 t>cc. Examples of the fits to the mean profiles are 
shown in Fig. 20. We note the following key points from 
this analysis. 

The overall quality of the line profile matches is 
good, which partly reflects the flexible optical depth ‘law 1 
adopted. It is notable that mostly in Civ - and to a 
lesser extent in Si IV - the low velocity absorption (i.e. 
below about 0.25ttoo) is not well matched. The spherically 
symmetric models often predict excess forward-scattered 
emission in the low-velocity part of the P Cygni absorp- 
tion trough. The discrepancy scales qualitatively with the 
strength of a given spectral line, such that it is gener- 
ally worse for the maximum absorption cases than the 
minimum ones for a single star. This does not hold as a 
function of mass-loss rate though, and for example the 
Si iv mismatch is just as large for HD 150168 as it is for 
HD 37128, for which M <?,- is a factor of ~ 5 larger (Fig. 20). 
Puis et al. (1993) discuss in detail the fact that reduced 
emission results when account is taken of back-scattering 
in several different resonance zones in structured winds. 
The reduced emission acts so as to yield an apparent in- 
creased absorption in the low velocity region between ~ 0 
and 0.3t>rx>. Also if the winds are not spherically symmet- 
ric, then the models will produce too much scattered light 
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. Table 4. Results from SHI model profile fits 


M q( N ) Mg(Si’ + ) A/r/(C"*+) XI q(Sl 2 + ) 


Star (HD) 

f'x 

mean 

range 

mean 

range 

mean 

range 

mean 

57 1 JS 

1700 

2.751-9) 

2.G0- 3. IO(-9) 

1.64 (-8) 

1. 13-l.79(-8) 

>4.79(-9) 

- 

l.56(-9) 

91969 

1500 

4.31 (-9) 

4.00 5.25(-9) 

1.60(-8) 

0.88-1.60(-8) 

>4.01(-9) 

- 

1.39(-9) 

164402 

1650 

6.31(10) 

5.2l-6.54t-10) 

lT2(-9) 

1. 12-1.66(-9) 

1.07(-9) 

0.82-1. 30(-9) 

1.34(-10) 

204172 

1750 

8.3G(- 10) 

0.71 1.53(-9) 

3.06(-9) 

2. 14-4. 34 (-9) 

1.92(-9) 

1.84-3.48(-9) 

3. 1 2 ( - 1 0 ) 

64760 

1500 

9.90(- 10) 

0.83- 1.131-9) 

2.88(-9) 

1.77-4.80(-9) 

1.82 (-9) 

1.54-2. 41 (-9) 

2.221-10) 

167756 

1750 

1.52(-9) 

1.40- l.C4(-9) 

3.25(-9) 

2.42-3. 37(-9) 

2.01 (-9) 

1.93—2.2 1 (-9 ) 

3. 70(- 10 ) 

150168 

1370 

5.30(-9) 

4.70-6.98(-9) 

3.58(-9) 

0.80-5. 76(-9) 

2. 34 (-9) 

1.55-2. 80 (-9) 

2 . 77 (-10) 

47240 

1000 

2.67 ( - 10 ) 

2.65 -3.87(-10) 

2.57(-9) 

2.20-5.43(-9) 

> 1 . 28( -9 ) 

- 

6.49(- 10) 

157246 

600 

7.30(-l0) 

7.00-14.30(-10) 

9.80(- 10) 

5.70-L6.50(- 1 0) 

>4.61(-10) 

- 

- 

96248 

800 

- 

- 

l.-16(-9) 

1.45-1. 46(-9) 

5.81 (-9) 

4.29-6.03(-9) 

- 

53138 

520 

- 

- 

2.78(-9) 

2.55-3.07(-9) 

7.27(-10) 

-i 

to 

O 

5 



The products of mass-loss rate and ionization fraction (M q t ) are derived from fits to the mean profiles. (X.XX(-Y) = X.XX x 
10 ' A/q yr 1 ). The range in this vaiue determined from cases of overall minimum and maximum absorption is also listed. 
The IUE SWP images corresponding to the cases of minimum and maximum absorption that were modelled are: HD 37128 
(8100, 8130), HD 91969 (57187. 57426), HD 164402 (54200, 54225), HD 204172 (38957, 48914), HD 64760 (53522, 53709), 
HD 167756 (29001. 29049), HD 150168 (47098, 47198), HD 47240 (48953, 48930), HD 157246 (54172, 54103), HD 96248 
(57069, 57204), HD 53138 (30153, 30264). We also obtained a few results from fits to Aim - HD 204172, 2.72(-9); HD 47240, 

2. 17(-9); HD 96248. 3.38(-8) and Cii - HD 53138, 6.84(-10). 


at low velocities. Note, however, that the discrepancy at 
low velocities in a fast rotator like HD 64760 is no greater 
than that for an apparent moderate rotator of the same 
spectral type, HD 167756 (Fig. 20). It is also possible that 
to some degree the mis-match at low velocities reflects 
clumping in a highly structured wind so that parts of the 
stellar disk are being obscured by very different amounts 
of material than others. The low velocity mis-match is 
generally much lower for N v compared to Si iv and Civ 
(Fig. 20), implying then that the region occupied by the 
N V emission may be more uniformly filled. 

The mean values of M m listed in Table 4 are consis- 
tent with the dependence on luminosity for B supergiants 
reported by Prinjafe Massa (1997), who presented prelim- 
inary results for a wider, time-averaged survey. Typically 
the range in M qi implied by the variability between min- 
imum and maximum absorption line profiles is a factor of 
~ 1.5. A clear exception is the factor of ^ 7 change noted 
in M q(Si 3 +) for HD 150168; this star showed strong evi- 
dence for ionization variability in the time-series data (e.g. 
Fig. 12 and see below). This quantity describes the change 
in structures integrated over a velocity range from 0.2 to 
0.9 i'co , and is therefore intended tc reflect variations in 
the global mass-loss rate. However, the localised variations 
(i.e., in a given velocity bin) can be several times larger. 
For example, at 0.7 (cf. temporal variance in Fig. 18), 
the product of M q(Si 3 +) may vary between a factor of 
4 to 15 for supergiants in the spectral range BO to BO. 5. 
These changes reflect the degree of siructure in the wind. 

Our sample of stars is too small to derive useful infor- 
mation on the temperature dependence of the relative ion 
ratios. However, at a given effective temperature the val- 
ues of q{Si 3 + /q(N 4 +), or q{Si 3 +)/q(C 3 +), can differ by up 
to a factor of 5. We find that there are no significant differ- 


ences (or trends) in an individual star between the relative 
ion fraction derived from the mean profiles, minimum ab- 
sorption profiles and maximum absorption profiles. Over 
the spectral range BO to B 1 (la and lb) the mean ion 
ratios are q{Si 3+ )/q( N 4+ ) — 2 . 5 , q(Si 3 +)/q(C 3+ ) ~ 1.7 
and q(Si 2 +)/q(Si 3 +) - 0.1. HD 150168 and HD 157246 
(7 Ara) exhibit the most exceptional mean ion ratios of 
q(N* + )/q(Si 3+ ). In HD 150168 the ratio is almost 3 times 
larger than for other stars of similar parameters. The dif- 
ference is linked to the ionization change (of unknown 
origin) that occurred during the run (see Sect. 3.2); the 
increase in the ion ratio points to a temporal ionization 
shift from a 'normal state" to a high state’. The ratio of 
q(N 4 +)/ q (Si 3 +) in 7 Ara is also several times larger than 
for other Bis. This result may reflect a nitrogen rich stellar 
surface due rotational mixing (in a rapid rotator). As ex- 
pected, the q{C 3+ )/q(Si 3+ ) ratio of the BC star HD 96428 
is abnormal, and the mean ratio is between 4 to 8 times 
larger than the other stars in our sample. 

The behaviour of the ionic ratios as a function of ve- 
locity is shown in Fig. 21 for our program stars (based 
on fits to the mean profiles of the time-series). The 
q(Si 2+ )/q(Si 3+ ) ratio does not depend on abundance and 
its observed trend indicates a decrease in ionization state 
of the wind as a function of velocity. The relative ion mix- 
ture of q(Si 3 +)/<?(A r4 +) seems to be more stable between 
~ 0.4 to 0.8 (or ~ 2 to 5 stellar radii for the adopted 
empirical velocity law). 

HD 37128 (e Ori) is the only star in our sample with 
a measured radio flux, which may be used to provide 
a reliable estimate of the mass-loss rate. Blomme et al. 
( 2002 ) have conducted a radio and sub-mm study of e 
Ori, and they derive a value of M = 1.9 x 10 “ 6 M yr -1 . 
The (mean) results in Table 4 therefore give ion frac- 
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N V - mean C IV - mean Si IV - mean Si III - mean 

37128 


91969 


164402 


204172 


64760 


167756 


150160 


47240 


-10 12 -10 12 

Normalised Velocity (v/v„) 


Fig. 20. SEI model line profile fits for N V, Civ, Si IV, and Si III lines in our program stars. Examples are shown for the cases 
of constructed mean profiles (see also Table 4). 

lions for this star of q(N*+) ~ 1.5 x 10“ 3 , q(Si 3+ ) ~ to B1 of q(N 4 +) ~2 x 10~ 3 ,g(St 3+ ) ~5 x 10” 3 , q{C 3+ ) 
9.0 x 10 -3 . q(C 3+ ) > 2.6 x 10~ 3 , and q(Si 2 ^) ~ 8.0 x ~ 3 x 10~ 3 , and q(Sr+) ~ 1 x 10 -3 . The line-fitting 
10 - * 4 , Alternatively, using the theoretical mass-loss recipe analysis suggests therefore that none of the spectral lines 
of Vink et al. (2001) to determine the mass-loss rates for modelled here represent the dominant stage of ionization 
our target stars, yields mean ionization fractions over BO for early B supergiants (over the UV line-formation re- 
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gion). This is particularly surprising for the case of Si ,,+ 
and C' 54 " in the earliest B stars studied, such as t Ori. As in 
O-type stars, the other likely candidates l'or the dominant 
ion stages are Si 4+ . C 4+ . and .\ 3+ . 

It may be that the low ionizat ion fractions derived here 
are due to the effects of clumping in the stellar winds 
(e.g. Massa et ah 2002, in preparation). If the winds are 
principally composed of optically thick clumps separated 
by transparent ‘voids’, then each constant velocity sur- 
face will only partially contribute to the emission line. 
Furthermore, the face of the star wo ild not be completely 
covered by optically thick material at any velocity. This 
means the absorption part of the profile will not be very 
deep, thus also giving the impression of an unsaturated 
line. As a result, when fitting the spectral line with a ho- 
mogeneous wind model, we derive i finite optical depth 
(and lower values of M q t ) , even though the optical depth 
of the wind might be much larger if i he wind material was 
homogeneously distributed. 

6. Discussion 

We have shown that the stellar wind activity in early B 
supergiants is rich and varied. Figures 2 to 17 illustrate 
the fact that a selection of temporal signatures can occur 
in an individual star, and between stars. Different struc- 
tures can also occur in the same star at different epochs, 
though 7 Ara (HD 157246) is the on y extreme example of 
this behaviour in our sample. The physical connection be- 
tween these perturbations remains unresolved. There are 
two general forms of variability characteristics that appear 
to co-exist in the UV data in some star's: 

1. Firstly there are broad absorption features which have 
a large spatial coherence (assuming monotonic veloc- 
ity laws) and are relatively short-lived. Examples of 
these can be seen in HDs 64760 and 150168 as bowed 
features, and HDs 91969, 164402 204172 and 47240 as 
blue-ward accelerating absorption enhancements. 

2. In most of the stars listed above the broad absorption 
features are sometimes ‘accompanied 1 by DACs, which 
are characteristically different in that they have longer 
temporal coherence. These feal ures are much nar- 
rower in velocity, accelerate blue-wards more slowly, 
and are quickly confined to the highest velocities near 
the short-ward profile edge. Their origin is enigmatic 
partly because the UV time-series are mostly not ex- 
tended enough to monitor the recurrence time-scales 
(and properties) of these features. It is therefore diffi- 
cult to establish (or rule out) their connection to the 
extended absorption features (or bowed structures). 

The first of the two phenomena are usually linked with 
the effects of CIRs (e.g. Owocki et al. 1995; Cranmer 
and Owocki 1996) passing through the iine-of-sight. 
Specifically, localised features or in homogeneities at the 
stellar surface are predicted to alter the radiative force, 
which leads to the formation of high-density low-speed 
gas streams. These streams are then rammed by t he faster 
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Fig. 21. The relative ion ratios (q,/qj) fire plotted as a func- 
tion of velocity, based on SEI model fits to the mean line pro- 
files of our program stars. 

moving underlying’ wind as the star rotates, leading to a 
spiral pattern. The enhanced absorption identified above 
arises from an extended layer of near-uniform velocity that 
marks the initial response of the unperturbed wind to the 
slower CIR ahead. 

A further possibility then is that the slowly migrating 
DACs (item 2 above) are also related to CIRs, but orig- 
inate in the radiative-acoustic kinks which are predicted 
to trail the CIRs in the hydrodynamical simulations of 
Cranmer & Owocki (1996). 

It is important to realise that both observational phe- 
nomena highlighted above are not always present in every 
OB star. In many O-type stars for instance, only the long- 
lived DAC features are present, i.e. without evidence for 
coexisting ‘bowed’ structures or periodic absorption mod- 
ulations. In this respect, we believe the model developed 
by Hamann et al. ( 2001 ) is more appropriate for a star 
with a distinctly structured wind such as HD 64760 (Fig. 
10) or £ Per (Kaper et al. 1997), than for £ Pup where 
the dominant (19.2-hr) variability pattern is due solely to 
blue-ward migrating DACs. 

The strong ionization changes seen in HDs 150168, 
91969, 64760 and 157246 (Table 3) are likely shock sig- 
natures, though we cannot easily assign them to the same 
physical origin. If the interface between slow and fast wind 
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streams in the ( IR model is associated with strong shocks, 
then we can expect, some ionization changes in the UV 
lines as the structure crosses the line-of-sight to the stel- 
lar disk. In cases where the variability phenomenon is seen 
at different velocities in different ions (Table 3). this be- 
haviour may also relate to the action of spatially confined 
shock regions in the wind, where the post-shock gas is 
denser, lower speed, and less highly ionized than the pre- 
shock material. Multan (1984' id Cranmer k Owocki 

(1996) discuss scenarios whici av potentially result in 

velocity jumps of several hundred km s _l . 

An important result from our study is the variety of 
wind activity patterns that are present in the U V time- 
series data of Bis. However, this diversity of behaviour 
does not necessarily indicate that there is a correspond- 
ing variety of fundamental processes responsible for gen- 
erating wind structure. The combination of stellar rota- 
tion rate and viewing (aspect) angle to the observer are 
among the most variable factors in our sample of BO to B1 
stars; more so than e.g. T e ff, log(L*), i’oo and wind density. 
Inclination angle and v rot could play an important role 
in determining the precise observational manifestation of 
azimuthally-extended large-scale wind structure that has 
its roots at the photosphere. 

An interesting signature of spiral-shaped structure in 
a hot star wind is the appearance of bowed feat ures. Since 
the winding of a spiral is determined by radial expan- 
sion and stellar rotation, the visibility of bowed features 
in UV time-series data is largely determined by the ratio 
of t'rot A oo- Though we only know the projected rotation 
velocity, and our statistical sample is this paper small, the 
ratio of v € sin(i)/voo is largest for HD 64740 (~ 0.14) and 
HD 157246 (7 Ara; 0.41 for 1995 data set). HD 64760 of 
course exhibits the clearest examples of phase bowing in 
our UV sample (Fig. 10), and indeed it is the best known 
case for this behaviour out of all the UV time-series data 
available for OB stars. Fast rotation is also undoubtedly 
the key to 7 Ara's extreme behaviour, i.e. the transtent 
two-component (polar and equatorial) wind. Interestingly, 
HD 150168 is the only other star in our study that shows 
bowed structures which are similar to those evident in 
HD 64760 (Fig. 12). The projected rotation velocity of 
this star is not exceptional and v e sin(i)/vco is ^ 0.09. 
This is possibly then an example then where either the 
inclination angle and/or the density contrast of gas in the 
interacting spiral streams favours the detection of phase 
bowing. 

We have discussed in this paper the best UV time- 
series data sets of early B supergiants available in any of 
the public archives. With the demise of IU E , these col- 
lections will not be improved for the forseeable future. 
One of our key objectives was to collect these homoge- 
neous data sets into a single comparative study. Over 
the next few years these results will be complemented 
by extended time-series studies of luminous early-type 
stars in the Magellanic Clouds using the Far Ultraviolet 
Spectroscopic Explorer satellite (FUSE: AA900 to 1200). 
We anticipate in particular more sensitive searches for ion- 


ization variability in the far-UV lines, and an examina- 
tion of the dependence of wind variability on metallicity. 
Additionally. ~2 metre class optical telescopes are cur- 
rently being used to monitor time-dependent, behaviour 
over rotational time-scales in Balmer and He lines of OB 
supergiants. A key goal is these studies is to exploit the 
recombination-dominated lines to determine the coher- 
ence of large-scale wind st-ri: ure 111 the deeper, iv-ar- 

star wind regions, lose to potential photospheric inho- 
mogeneities. 
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